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1.0  INTRODUCTION 


The  many  ways  In  which  the  earth's  upper  atmosphere  and  Ionosphere  Influence 
military  hardware  systems  have  been  discussed  In  past  MDAC  reports  and  proposals 


militant  h 
to  onr/t 


to  ONR.  This  annual  report  is  concerned  with  work  In  progress  on  the  development 
of  a quantitative  global  model  of  the  density  of  electrons,  N^,  In  the  upper 
atmosphere.  Knowledge  of  the  ionosphere  electron  density  Is  essential  for  the 
quantitative  determination  of  radio  wave  propagation  through  this  region  of  th< 
earth's  environment. 


Ultimately  N^  will  be  described  on  a completely  physical  basis.  This  approach 
to  the  electron  density  problem  has  been  taken  by  NRL  and  should  eventually 
provide  the  best  models  of  n£  The  purpose  of  the  current  contract  has  been  to 


provide  the^best  model 


} -currently-  for  the  electron 


density.  In  order  that  this  model  be  useful  to  groups  whose  function  is  impacted 
by  Ionospheric  variability  on  comnunlcatlons,  the  model  must  meet  several  require- 
ments. First,  It  must  be  global  (there  have  been  several  models  developed  of 
for  a given  latitude  or  local  time  region.  They  all  encounter  difficulties  when 
attempts  are  made  to  link  them  with  models  describing  other  regions).  The  mode; 
must  be  analytic*  There  are  two  reasons  for  this.  First,  an  analytic  model  will 
almost  certainly  be  capable  of  being  described  in  a computer  code  that  is  fast 
running.  Second,  the  model  almost  has  to  be  analytic  If  It  Is  to  be  used  to 

describe  In  real  time  changes  In  n£.  It  Is  also  prudent  to  use  an  analytic 

4 

model  to  determine  total  electron  content  (TEC)  by  Integrating  Nfi  along  some 
path.  The  model  Is  Intended  only  to  describe  gross  features  of  the  Ionosphere 


*v  ■%* v:  *rt-  v 


and  their  tempera!  variations.  It  is  not  capable  of  describing  the  more  isolated 

* 

rapid  fluctuations  In  N'  that  occur  in  some  regions  of  the  ionosphere.  * 

The  Intended  form  of  the  mode!  may  be  further  described  by  comparing  It  with 
other  Ionospheric  models.  The  NRL  one-dimensional  model  Is  meant  to  be  based  on 
a rather  complete  understanding  of  the  physics  and  chemistry  of  the  upper  atmo- 
sphere and  Ionosphere.  The  model  determines  atmospheric  density  and  winds  and 
uses  a detailed  knowledge  of  chemistry  of  the  upper  atmosphere  together  with 
other  physical  phenomena  to  finally  determine  at  a given  latitude  the  variation 
In  Ng  and  local  time.  This  model  takes  a fairly  large  amount  of  computer  time 
and  cannot  be  considered  usable  for  real  time  specification  and  prediction  of 
Nfi.  However  It  Is  clearly  one  step  In  the  approach  toward  the  ultimate  3D 
physical  model  (which  also,  by  virtue  of  Its  complexity,  may  not  be  useful  for 
the  solution  of  real  time  problems).  The  Aerospace  Corporation  has  developed 
an  Ionospheric  model  that  is  analytic  but  Is  merely  a functional  representation 
of  Nfi.  In  Its  construction,  little  consideration  was  given  to  the  relation 
between  the  quantitative  model  and  the  Input  parameters  it  might  use.  This 
model  Is  therefore  quite  limited  In  the  sense  that  It  cannot  be  used  In  a 
predictive  mode.  Its  associated  computer  code,  however,  is  much  faster  than 
the  NRL  model.  Work  begun  at  Air  Force  Global  Weather  Central  and  continuing 
now  at  the  Air  Force  Geophysics  Laboratory  describes  a statistical  Ionospheric 
model.  Such  a model  will  be  of  use  to  GWC  In  making  average  dally  forecasts 
but  will  be  of  limited  value  for  the  real  time  specification  of  prediction  of 
Nfi.  In  the  modeling  exercise  described  here,  an  attempt  has  been  made  to  meet 
all  of  the  requirements  discussed  above  (that  the  model  be  global,  fast,  analytic. 


2 


■v. 


Wt 


r 


i 


and  describe  all  gross  global  Ionospheric  features).  This  has  dictated  that  the 
model  be  unlike  all  of  the  others  just  discussed.  Instead,  it  is  semlemplrlcal , 
based  on  as  much  physics  as  can  be  used  to  accurately  and  quickly  describe 
features  of  Nfi.  Also,  a great  deal  of  attention  has  been  paid  to  the  relation 
between  the  model  and  the  sources  of  data  that  might  be  Input  to  It.  Any  model 
Is  really  no  more  accurate  than  the  Input  data  supplied  to  It.  These  inputs 
hopefully  will  Include  Information  provided  from  sensors  on  the  SOLRAD-HI, 

DMSP,  and  HELIOCENTRIC  satellites.  Such  data,  we  believe,  will  result  In  a much 
more  accurate  description  of  Ng  than  can  be  provided  merely  by  use  of  the  10.7  cm 
flux  Index  and  magnetic  Indices  such  as  K . This  permits  the  model  to  use  real 
time  data  and  not  Indices  that  are  averaged  over  space  and  sometimes  not  available 
for  at  least  minutes  and  sometimes  days  and  weeks  after  the  time  of  interest. 

Some  of  the  data,  for  example  the  Interplanetary  and  solar  wind  Information 
provided  by  the  HELIOCENTRIC  satellite,  may  permit  the  prediction  of  some  of 
those  high  latitude  Ionospheric  features  under  the  control  of  low  energy  charged 
particles  emanating  from  the  sun.  This  might  provide  "lead  times"  of  up  to  one 
half  hour.  Other  changes  In  the  density  of  the  neutral  atmosphere  and  in  Ng  at  low 
latitudes  caused  by  events  initiated  In  the  polar  and  auroral  regions  may  be 
predictable  with  lead  times  of  several  hours. 

In  the  model  both  "primary"  and  "secondary"  parameters  are  used.  The  primary 
parameters  are  those  whose  temporal  variations  are  described  directly  by  Input 
data  to  the  model.  The  set  of  secondary  parameters  are  constants  used  to 
describe  several  Ionospheric  features.  These  constants  do  not  vary  with  the 
Input  parameters.  Rather  they  have  been  determined  from  various  observational 
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data  sets  describing  particular  gross  ionospheric  features.  A particular 
feature  In  the  model  then  Is  first  described  in  some  average  sense  In  terms 

t 

of  secondary  parameters.  Temperal  and  spatial  variability  of  a particular 
feature  is  determined  by  the  changes  in  its  primary  parameters  which  are  In 
turn  controlled  directly  by  input  data  to  the  model. 


Not  all  of  the  Input  parameters  will  link  the  model  directly  to  spacecraft 
data.  Instead,  some  of  them  will  be  tied  to  other  MDAC  models  of  environmental 
features.  Most  Important,  the  model  of  Ng  relies  directly  on  the  MDAC  model 
of  neutral  atmospheric  density.  Conceptually,  the  electron  density  may  be 
thought  of  as  a product  of  the  neutral  atmospheric  density  times  some  flux 
function.  Actually,  the  neutral  density  model  Is  used  to  integrate  the  flux 
of  electromagnetic  energy  from  the  sun  down  through  the  atmosphere  in  order  to 
determine  Ionization  profiles.  Use  Is  also  being  made  of  the  MDAC  magneto- 
spheric  magnetic  field  model.  This  is  because  many  Ionospheric  features  are 
under  control  of  the  geomagnetic  field.  In  particular,  several  of  the  high 
latitude  dayslde  features  vary  as  the  dayside  cleft  (or  dayslde  cusp)  regions 
change  their  position.  These  other  MDAC  models  themselves  have  Inputs  from 
satellite  data.  For  example,  the  neutral  density  model  relies  on  data 
describing  the  solar  electromagnetic  spectrum  and  other  parameters  describing 
charged  particle  energy  sources  at  high  latitudes.  A detailed  description  of 
the  model  Is  provided  In  the  following  section  where  the  geomagnetic  control 
of  some  Ionospheric  features  and  this  model's  relations  to  other  MDAC  models 
are  discussed. 


As  mentioned  above  It  Is  our  Intent  that  this  model  provide  the  current  most 
accurate  description  of  Ng  In  some  global  sense  and  that  tie  model  be  con- 
structed such  that  It  Is  fast  on  the  computer  and  capable  of  describing  Ng  In 
both  the  specification  and  prediction  modes.  We  anticipate  that  the  model  will 
be  used  Immediately  for  both  systems  applications  and  scientific  purposes.  We 
would  hope  that  It  might  be  used  with  our  neutral  density  model  to  provide  the 
first  accurate  quantitative  descriptions  of  global  atmospheric  wind  systems 
(a  major  difficulty  concerning  workers  on  atmospheric  winds  Is  the  lack  of  a 
good  description  of  both  neutral  density  and  electron  density  throughout 
the  upper  atmosphere).  It  Is  also  our  opinion  that  the  development  of  the 
model  has  provided  several  Insights  into  the  linking  of  the  Ionosphere  with 
the  upper  atmosphere  and  magnetosphere  and  that  the  availability  of  such  a 
model  should  hasten  the  development  of  a more  sophisticated  second  generation 
quantitative  description  of  ionospheric  electron  density.  Through  the  develop- 
ment of  this  model  we  have  been  In  contact  with  NRL,  NELC , AFGWC,  AFWL,  and 
Aerospace  Corporation.  These  groups  will  be  among  the  first  users  of  the  model. 

We  hope  to  have  a large  amount  of  feedback  from  them.  We  have  found  with  our 
other  models,  especially  the  MDAC  magnetospherlc  magnetic  field  model  (developed 
with  ONR  funding),  that  this  Interaction  with  model  users  Is  vital.  It  has  In 
the  case  of  the  magnetic  field  model  led  to  a much  more  accurate  second  generation 
model  that  has  just  recently  been  completed. 

We  are  proceeding  as  planned.  The  program  plan  for  the  past  year  called  for 
determining  Nfi  for  the  several  global  Ionospheric  features  we  planned  to  represent, 
combine  those  features  to  form  the  global  model,  and  choose  provisional  Inputs  fir 


5 


■,  »•»»  •• . 


v.Hk. 


*v  ‘M'./WV  ' 


l 


the  model.  These  three  tasks  have  been  completed  and  the  resulting  model  is 


described  in  the  next  section.  The  Important  task  of  "model  calibration" 
remains  and  it  is  planned  that  this  work  will  be  Initiated  in  1977.  It  is  this 
calibration  that  will  tie  the  model  to  the  satellite  data  sets  and  make  it,  we 
hope,  a very  useful  tool  for  specifying  and  predicting  global  features  of  Ng. 

Details  of  the  form  of  the  model  are  presented  in  the  following  section.  We 
plan  to  present  a paper  on  this  model  at  the  Spring  American  Geophysical  Union 
Meeting  and  at  other  meetings  as  appropriate.  We  further  plan,  once  the  cali- 
bration is  completed,  to  test  the  model  against  several  observed  structures  in 
Ng  and  finally  to  publish  papers  describing  its  construction,  form,  and  uses. 


2.0  MODEL  DESCRIPTION 


We  report  here  the  construction  of  the  model  In  terms  of  Its  primary  and  secondary 
constants  and  parameters  (see  Section  1).  The  model  specifically  includes:  the 

structure  of  the  E,  Fj , and  F2  layers,  the  equatorial  anomaly,  the  winter  anomaly, 
polar  cap  structure,  semiannual  and  seasonal  dependences,  and  the  mid-latitude 
trough. 


The  input  parameters  to  the  model  are  chosen  such  that  satellite  data  can  be 
used  together  with  the  model  to  specify  Ng  in  real  time.  Hopefully,  with  data 
from  the  HELIOCENTRIC  satellite,  the  model  can  be  used  to  predict  those  features 
controlled  by  magnetic  fields  and  solar  charged  particles.  In  part,  the  pre- 
dictive capability  of  the  model  stems  from  its  dependence  on  Inputs  from  other 
quantitative  environmental  models.  For  example,  the  location  of  the  mid- latitude 
trough  Is  largely  controlled  by  the  locations  (L  values)  of  the  plasmapause  and 
the  auroral  oval.  Coupling  of  the  model  to  the  MDAC  magnetic  field  model  pro- 
vides the  necessary  Input  on  the  dayslde  cusp  position.  More  important,  the 
model  relies  heavily  on  the  MDAC  model  of  upper  atmospheric  neutral  density,  N. 
Several  Ionospheric  structures  and  temporal  dependences  are  found  to  tie  direct  y 
to  variations  in  atmospheric  density.  This  is  Illustrated  In  Figure  1 where  the 

semiannual  variation  In  N,  at  the  F,  peak  is  shown.  This  variation  In  N,  is 

e c.  e 

produced  entirely  by  the  semiannual  variation  In  N described  In  the  MDAC  neutral 
density  model.  In  addition,  the  structures  of  the  E and  F1  layers  are  determined 
by  calculating  the  absorption  of  solar  electromagnetic  radiation  In  the  neutral 
density  model. 


...  ■.  ♦ .. ...  . t ^ . ... . ,■»  - 


^ ^ .'*K.***.v 


V : * .* 


ANNUAL  VARIATION  AT  THE  F2  PEAK 


DAY  OF  YEAR 


The  semi  empirical  formulation  of  the  model  is  discussed  below  with  emphasis  on 
layer  formation  and  anomalies. 
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2.1  Layer  Formation 

Let  S be  the  energy  flux  at  location  I and  S + dS  the  flux  at  location  i + dt. 

Let  o be  the  absorption  cross  section  of  the  atoms  of  the  gas  and  N their  number 
density  (the  neutral  density  of  the  upper  atmosphere).  The  energy  absorbed,  dS, 
in  a cylinder  of  unit  cross  section  and  axis  parallel  to  the  direction  of  the 
incident  beam  is  given  by 

dS  = SoN  df. 

• ■; 

Nd£  = - t 

where  x is  the  optical  thickness  of  the  atmosphere  at  location  ZQ  along  the  path 
to  the  sun.  As  S S (its  value  above  the  atmosphere)  x -*■  0. 

oo 

Thus  S = Sm  exp  (-x), 

and  the  energy  absorbed  per  unit  volume  is 

gf-  ■ NoS  = HoSm  exp  (-x). 

If  n is  the  production  efficiency  per  unit  energy,  then  q (the  number  if  ion/ 
electron  pairs  produced  per  unit  volume)  Is  given  by 

qU)  ■ NanS^  exp  (-x). 
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The  above  equation  gives  the  production  of  electrons.  The  two  predominant  loss 
mechanisms  for  lon/electron  pairs  are  recombination  and  attachment.  For 
recombination  then 


e + 

dr  * Ne  N 


where  is  the  electron  density  N the  positive  ion  density  and  a the 


recombination  coefficient. 


Thus,  at  equilibrium 


ar  ■ >>  Ne N*  ■ 0 


Then,  assuming  that  N £ N , 


Non  S exp  ( -t ) -a  N d = 0 


Non  1/2 

Therefore,  Ng  * [ — exp  (-t)] 


For  attachment  the  loss  equation  is 


dN 

r=-bNeN 


» *♦  ***  + **  • ^ ** ' 
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In  discussing  attachment  processes  In  the  Ionosphere  it  is  generally  assumed 
that  the  number  of  neutrals  Is  much  greater  than  the  ions  so  that  when  Ionization 
occurs  the  neutral  density  does  not  change  significantly.  Thus 


ar 3 Ne 


Is  generally  used,  where  3 is  the  attachment  coefficient. 


Then,  at  equilibrium 


= q -6  Nfi  = 0 and 


Non  „ 

= — g—  exp  (-t) 


In  general  a and  6 may  vary  with  height  because  the  reactions  usually  involve 
three  bodies  Instead  of  two. 

2.1.1  The  E and  F1  Layers 

It  has  been  observed  that  the  E and  the  F-j  layers  of  the  ionosphere  can  be 
explained  quite  well  using  the  recombination  equation.  Thus,  for  the  E and  F.j  ayers 


Ng  3 [N  SE  aE  exp  ( -te ) ] 

Fl  1/2 

N - [N  SF  dp  exp  (-Tp  )] 
e F1  F1  F1 


*■  - ■ « i V»  1 , . »^. ..  ...  • 


E E 

where  N and  N are  the  electron  concentrations  for  the  E and  F,  layers,  S 
e e i 

and  S 1 are  the  solar  U.V.  flux  affecting  the  E and  layers,  a£  and  are 
the  recombination  coefficients,  N the  neutral  density,  and  given  by 


N dt  £ ot 


N dl 


where  T = 


N d£  is  the  total  atmospheric  cross  section  between  the  observation 


point  and  the  sun.  The  MDAC  atmospheric  density  model  (which  Is  analytic  and 
computationally  very  fast)  is  used  in  the  evaluation  of  the  above  integral. 


Then  TE  = aE  ^ 
and  = ar  T 

n hl 

where  or  and  oc  are  the  absorption  coefficients  for  the  wavelengths  affecting 
h hl 

the  E and  F1  layers. 

These  equations  describe  the  observed  structures  in  the  E and  F^  layers.  Chang »s 
In  T depend  on  the  neutral  density  and  the  location  of  the  sun.  Thus,  by  using 
the  MDAC  neutral  density  model  (which  contains  most  of  the  variations  observed 
In  the  neutral  atmosphere)  and  a simple  sun-position  program,  the  variations 
observed  In  the  E and  F^  layers  can  be  accurately  described. 

The  Fj  layer  Is  essentially  turned  off  at  local  sunset  but  the  E layer  is  maintained 
through  the  night.  Thus 


r ' 


r . . 1/2 

Ng  (night)  “ [N  S£  a£  exp  (-te)] 

I 

where  SE  Is  the  Intensity  of  scattered  light  and 
*o 

te  * aE  | N di 
00 

where  j N di  Is  the  optical  thickness  (Integrated  vertically  - not  toward  the 
sun). 


Contours  of  Ng  at  noon  and  midnight  are  shown  at  the  equator  and  mid  latitudes 
In  Figure  2.  These  results  are  obtained  simply  by  Integrating  the  solar  flux 
through  the  neutral  density  model.  The  model  Is  constructed  so  that  as  param- 
eters describing  the  solar  flux  vary  in  time  (for  example,  as  observed  by 
SOLRAD-HI)  the  Input  parameters  will  vary  such  that  the  output,  Nfi,  will  exhibit 
the  proper  temporal  variation. 


2.1.2  The  Fq  Layer 

The  structure  and  variability  in  the  F2  layer  is  not  so  directly  controlled  by 
the  neutral  atmosphere.  All  attempts  to  treat  the  F2  layer  using  U.V.  generate^ 
Ionization  and  electron  losses  which  are  a function  only  of  density  ended  In 
failure.  It  Is  generally  agreed  that  the  bottom  of  the  F2  layer  Is  defined  as 
the  region  where  electron  loss  rates  change  dramatically  (from  the  order  of 
seconds  to  the  order  of  hours).  Several  studies  were  undertaken  with  neutral 
constituents  and  electron  density  dependent  attachment  and  recombination  coeffi- 
cients. In  none  of  these  attempts  was  It  possible  to  reproduce  the  observed 
seasonal,  semi-annual  or  solar  cycle  dependencies. 
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ELECTRON  CONCENTRATION  (cm 


An  empirical  function  was  therefore  developed  which  arbitrarily  (adjusted  to 
provide  a best  fit  at  some  point  In  time)  defines  a loss  term  varying  with 
height  and  an  Ionization  term  depending  only  on  the  ambient  density.  Further 
empirical  observations  suggested  that  the  normal  layer  was  limited  to  the 
region  equatorward  of  the  latitude  of  the  plasmasphere.  Thus,  a function  was 
Introduced  which  limits  the  F2  layer  (except  for  the  winter  anomaly  peak)  to 
L values  ^ 4-6.  (The  L value  can  roughly  be  thought  of  as  the  equatorial 
extent  of  magnetic  field  lines  measured  from  the  center  of  the  earth  In  earth 
radii.  Thus  at  the  equator  L = 1.)  That  Is,  the  neutral  density  dependent 
F2  electron  values  are  multiplied  by  the  function 


where  A * cos 


L( A)  * exp  (-y  AJ) 
1 


(V 

V L . 


Is  the  magnetic  Invariant  latitude.  The  parameter  y 
was  adjusted  to  give  a best  fit  to  the  midnight  equatorward  edge  of  the  mid 
latitude  trough.  Thus,  the  parameter  y is  a primary  parameter,  since  it  is 
derived  from  the  size  of  the  plasmasphere.  To  properly  calibrate  the  parameter 
Y simultaneous  satellite  measurements  of  the  plasmapause  and  location  of  the 
equatorward  edge  of  the  trough  are  required. 


The  Fg  layer  Is  essentially  Independent  of  direct  solar  zenith  angle  control, 
and  depends  Instead  on  the  density  of  the  neutral  atmosphere  and  the  latitude 
extent  of  the  plasmasphere.  The  winter  anomaly  however  is  observed  to  be 
under  direct  solar  control.  It  appears  at  sunrise  and  rapidly  disappears  at 
sunset.  Much  of  the  average  structure  of  the  F2  region  as  given  by  the  model, 
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Including  the  peak,  is  shown  In  Figures  3 and  4 In  terms  of  the  critical 
frequency  (which  relates  directly  to  Ng).  Note  that  on  the  noon  meridian 
some  of  the  structure  Is  produced  by  the  equatorial  anomaly  (around  the 
12  Mhz  contour).  The  magnetic  latitude  dependence  of  the  equatorial  anomaly 
Is  Included  In  the  model  and  discussed  below. 


2.2  Equatorial  Anomaly 

The  equatorial  anomaly  which  is  thought  to  be  caused  by  an  upwelllng  in  the 
equatorial  region  produced  by  electric  fields  may  be  organized  around  a 
central  L value,  LQ.  The  functional  form  of  the  anomaly  term  is 


d-exp  C-5(L-L0) ] 


where  the  amplitude,  d,  is  dependent  on  solar  cycle  as  well  as  local  time  and 
6,  which  defines  the  width  of  the  anomaly  peaks,  has  a local  time  and  solar 
cycle  dependence  and  a latitude  dependence.  The  anomaly  Is  strongest  along 
the  magnetic  field  line  with  the  L value,  LQ.  The  anomaly  dependence  on 
altitude  and  magnetic  latitude  as  given  by  the  model  are  shown  In  Figure  5. 


2.3  Winter  Anomaly 

In  this  model  the  winter  anomaly  is  treated  as  a completely  separate  source  and 
loss  mechanism  and  Its  contribution  to  the  electron  content  Is  additive  to  all 
other  electron  sources.  The  formalism  developed  for  the  F^  layer  was  used  for 
electrons  associated  with  the  winter  anomaly.  The  winter  anomaly's  contribution 
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2W 

to  the  F2  layer,  Ng  , 
written  as 


during  the  winter  months  near  solar  maximum  is  therefore 


F2W  1/2 

N - g(A,P)  [NSp  Or  exp(-Tp  )] 

e r2W  2W  r2W 


where  g(A,p)  Is  a latitude  and  solar  cycle  dependent  function  which  determines 

the  extent  and  amplitude  of  the  winter  anomaly.  N is  the  neutral  density,  SP 

Is  the  solar  flux  responsible  for  the  Ionization,  aP  Is  the  Ionization 

r2W 

efficiency,  and  tf  ■ oP  T,  where  oF  is  the  absorption  coefficient  and  T 
r2W  r2W  r2W 

the  optical  density  to  the  sun. 


In  general  there  Is  no  F2  layer  over  the  polar  cap  In  the  above  formulation 
(except  for  the  winter  anomaly  which  at  times  extends  into  the  polar  cap).  The 
Fj  layer,  however,  extends  to  higher  and  higher  altitudes  as  the  light-dark 
terminator  Is  reached. 

j 

■ 

2.4  Polar  Cap  Structure  and  Mid-Latitude  Trough 

To  complete  the  polar  cap  structure  the  charged  particle  source  must  be  included. 
These  high  latitude  sources  appear  most  often  at  auroral  latitudes  (at  L values 
* greater  than  8.5-9).  At  night  the  particles  originate  from  the  tall  and  are 

more  energetic  than  the  particles  which  enter  through  the  dayside  cusp  regions. 

j 


The  configurations  of  the  auroral  E and  F layers  vary  extensively  In  amplitude 
and  location.  For  sake  of  simplicity  a simple  Gaussian  has  been  used  centered 


» 

* 


at  L a 10.5  with  a low  latitude  edge  of  8.5  at  midnight.  During  daylight  hours 
the  low  latitude  edge  moves  to  slightly  higher  latitudes  and  Is  centered  on  the 
dayslde  cusp  field  line. 


The  disappearance  of  the  F2  plasmasphere  associated  electron  density  together 
with  the  corpuscular  generated  auroral  electron  layer  forms  the  midlatitude 
trough.  During  daylight  hours,  especially  during  the  winter,  the  U.V.  produced 
electron  concentration  may  "fill  in"  the  trough.  The  variation  In  the  F2  peak 
through  the  trough  region  Is  shown  for  the  "average"  model  In  Figure  6.  Note 
that  the  important  feature  In  Ng  Is  largely  under  magnetic  control.  That  Is, 
as  the  plasmasphere  expands  or  contracts,  the  low  latitude  edge  of  the  trough 
moves  respectively  to  higher  or  lower  latitudes.  Likewise  the  field  lines 
along  which  the  dayslde  cusp  particles  precipitate  can  also  change  their 
latitude  and  simultaneously  the  latitude  of  the  poleward  edge  of  the  trough. 


ELECTRONS  (107cm' 


3.0  COMMENTS  AND  SUMMARY 


In  the  above  description  little  was  said  concerning  the  primary  and  secondary 
parameters.  This  Is  primarily  because  the  parameters  will  be  determined  by 
satellite  and  other  real  time  data  sets.  The  values  of  these  parameters  will 
be  listed  later  when  an  "average"  model  Is  published.  For  example,  the  value  of  LQ 
(used  to  organize  the  equatorial  anomaly  information)  currently  Is  set  at  1.12. 
During  the  next  year,  when  the  model  calibration  Is  being  performed,  this  value 
may  be  changed.  More  Important,  the  primary  parameters  associated  with  the 
equatorial  anomaly  determination  will  be  determined.  They  will  provide  the 
model  with  Its  predictive  capability. 


It  Is  emphasized  again  that  the  model  is  not  statistical  and  that  when  the 
calibration  Is  completed.  It  will  be  possible  to  use  It  in  real  time  to  specify 
and  in  some  cases  predict  temporal  and  spatial  variations  in  Ng.  It  Is  intended 

I 

that  this  model  be  used  with  satellite  data  sets  (although  some  provision  may  bo 
made  for  Its  use  with  magnetic  Indices)  and  other  quantitative  models  of  the 
near  earth  environment  such  as  those  developed  by  the  MDAC  group.  Further  it 
Is  our  hope  that  this  model  will  be  of  use  to  both  systems  groups  and  research 

» 

people.  The  Immediate  model  uses  with  communications  problems  are  obvious.  It 
Is  hoped  that  the  model  will  also  be  used  for  such  research  endeavors  as  the  deter- 
mination and  understanding  of  the  global  upper  atmospheric  wind  system. 

I 
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